Temperature dependence and anisotropy of the bulk upper critical field H c2 of MgB 2 
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The bulk upper critical field H C 2(T) of superconducting MgB2 and its anisotropy are established 
by analyzing experimental data on the temperature and magnetic field dependences of the afe-plane 
thermal conductivity of a single-crystalline sample in external magnetic fields oriented both parallel 
{Hc2) and perpendicular {H& ) to the c axis of the hexagonal lattice. From numerical fits we deduce 
the anisotropy ratio 70 = (0) /H^ 2 (0) = 4.2 at T — K. Both the values and the temperature 
dependences of H^ 2 and H^2 are distinctly different from previous claims based on measurements 
of the electrical resistivity. 

PACS numbers: 74.60.Ec, 74.70.-b, 74.25.Fy 
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Since the recent discovery of superconductivity of 
MgB2 with a critical temperature T c ~ 40 K, 1 a large 
number of experimental results on different properties of 
this compound has been reported in the literature. Most 
experiments were made using powder or polycrystalline 
samples. The hexagonal crystal structure of MgB 2 , how- 
ever, is expected to cause pronounced anisotropies in the 
electrical and magnetic properties, which can unambigu- 
ously be probed only by experiments using single crystals. 
In particular, the upper critical field H C 2(T) is an impor- 
tant parameter for characterizing the superconducting 
state of type II superconductors. 2,3 ' 4 For anisotropic ma- 
terials, such as hexagonal MgB2, the values of H C 2 may 
vary considerably for different orientations of the exter- 
nal magnetic field H. Choosing the field directions either 
perpendicular or parallel to the c-axis, the anisotropy 
may be expressed by a parameter 7 — H^ 2 /H^ 2 , which, 
in the most general case, may be temperature dependent. 
Earlier experimental results, mainly based on measure- 
ments of the electrical resistivity p(T,H), have resulted 
in a broad range of values of 7 and of extrapolated zero- 
temperature values of H^ 2 (0) and H® 2 (0) (for a review, 
see Ref. 5). Most of these experiments, also on single 
crystals, 6 ' 7,8 indicate a positive curvature of H c2 (T) in a 
wide range of temperature below T c and correspondingly, 
rather high critical fields at T = 0. Attempts to explain 
these features have lead to theoretical work suggesting 
the existence of some soft bosonic modes 2 and even un- 
conventional mechanisms of superconductivity have been 
considered. 9 

In this paper we present an evaluation of H^ 2 (T) and 
H$(T) of single crystalline MgB 2 , based on measure- 
ments of the thermal conductivity n(H,T). Comple- 
mentary results of p{T, H), obtained on the same single- 
crystalline sample, indicate that electrical transport mea- 
surements are not well suited to probe the bulk upper 
critical field H C 2(T) of MgB 2 . Inspecting the tempera- 
ture dependences of both H^ 2 (T) and H" 2 (T) close to 
to T c , our results indicate that even the zero-field crit- 
ical temperature T c (0) of the bulk may be lower than 
commonly believed up to now. This indicates that, in 
relation with superconductivity of MgB 2 , surface effects 
must be considered. 
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FIG. 1: The low temperature electrical resistivity of MgB2 
and its magnetic field dependence for a current in the ab- 
plane. The broken lines for H — 50 kOe indicate how the 
onset of the resistive transition defining Hp has been estab- 
lished. 



The thermal conductivity was measured in the basal 
plane of hexagonal MgB 2 exposed to varying magnetic 
fields H, oriented parallel and perpendicular to the basal 
afe-plane with a small misalignments of 3.5 ± 0.5° be- 
tween the field directions and the orientation of the plane. 
A standard uniaxial heat flow method, as described in 
Ref. 10, was used for the k(H, T) measurements. The 
temperature difference between the two thermometers 
was about 1% of the absolute average temperature. The 
measurements of the electrical resistivity p(H, T) were 
made using a 4-contact scheme and a dc-current of den- 
sity 50 A/cm 2 in the a6-plane with H along the o 
direction. The investigated single crystal has lateral di- 
mensions of 0.5 x 0.17 x 0.035 mm 3 and was grown em- 
ploying a high-pressure cubic anvil technique as described 
elsewhere. 11 

Low-temperature p(T) curves measured in constant 
external magnetic fields H are presented in Fig. 1 for 
T < 50 K. The zero- field resistive superconducting tran- 
sition at T c = 38.1 K is rather narrow (AT C ~ 0.15 K), 
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but the application of magnetic fields broadens the tran- 
sition considerably. If the field dependence of the on- 
set of the resistive transition, as illustrated in Fig. 1 for 
H = 5 T, is plotted in an [H,T] diagram, the curve de- 
noted as Hp in Fig. 4 is obtained. These H c p data are 
qualitatively and quantitatively very similar to results 
previously obtained on single crystals, 6 ' 7,8 in particular 
with respect to the positive curvature of H c p (T). In these 
earlier works H p (T) was associated with H^ 2 (T). 

The if-dependence of the thermal conductivity was 
measured at selected constant temperatures in the range 
between 2 and 50 K and in fields up to 60 kOe. Repre- 
sentative k(H) curves at selected temperatures are dis- 
played in Figs. 2 and 3 for H \\ c and H |j ab, respectively. 
As demonstrated in the inset of Fig. 2, a hysteretic be- 
havior of k(H), caused by vortex pinning, is observed in 
the low-field regime for H \\ c. In order to avoid ambigui- 
ties, each new field setting at a constant temperature was 
achieved by heating the sample to the normal state above 
50 K, and subsequently cooling it to the set temperature 
in the chosen field. In this way, a smooth variation of 
k(H), as demonstrated by the open circles (FC) in the 
inset of Fig. 2, was obtained. The field values H? r , be- 
low which the irreversibility is discernible, are rather low. 
For T = 4.03 K, e.g., if£ r - 0.7 kOe. At elevated tem- 
peratures and for H \\ ab the irreversibilities are reduced, 
as demonstrated in the inset of Fig. 3. 

The curves presented in Figs. 2 and 3 reveal the gen- 
eral features observed at all temperatures below 38.1 K. 
Starting at H = 0, n drops with a steep slope and, af- 
ter passing through a minimum, increases again until a 
region of very weak field dependence above some critical 
field, denoted as H K is reached. It is remarkable that, 
for each temperature, H^ is distinctly lower than Hp and 
that no distinct feature of k(H) is observed in the re- 
gion of Hp. This is explicitly demonstrated in Fig. 2. 
With increasing temperature, H K decreases towards zero 
as T approaches T c . This general n(H) features are typ- 
ical for type II superconductors and can be explained 
as follows. 12 The thermal conductivity of a superconduc- 
tor is due to itinerant electrons (n e ) and phonons (re p h)- 
Enhancing H from zero eventually causes the formation 
of vortices in the bulk of a type II superconductor. Af- 
ter zero-field cooling, the first vortices form at the lower 
critical field H c \. Consequentially, some additional scat- 
tering of phonons by normal electrons in the cores of the 
vortices will reduce K p h- With further increasing field the 
decrease of K p h is compensated by an enhancement of n e . 
Above H c2 , in the normal state, the field dependence of 
both K p h and n e is expected to be weak. The overall be- 
havior of the k(H) curves shown in Figs. 2 and 3 reflects 
these expectations, and as may be seen, k(H) is virtually 
field independent for H > H K . 

A more complete analysis of the k(H) data will be pre- 
sented in a forthcoming paper. 13 Here, we concentrate on 
the opportunity that these data allow for a reliable eval- 
uation of the bulk upper critical field H c2 (T), which ob- 
viously coincides with H K (T) as derived from our k(H) 
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FIG. 2: The magnetic field dependence of the thermal con- 
ductivity k(H) for H || c at T = 4.03, 7.93 and 20.5 K. The 
solid vertical arrows mark H%, set equal to the upper critical 
field Hc 2 - The dotted vertical arrow denotes Hp (see Fig. 1). 
The inset demonstrates the irreversible behavior of k(H) be- 
low Hf rl for T =4.03 K. 



curves for both field orientations. It may be seen that 
H^ 2 (T) = H%{T) is distinctly different from H c p (T). The 
solid line in Fig. 4, representing a general prediction for 
H C 2{T) of a conventional type II superconductor in the 
case where the coherence length £ and the electron mean 
free path I are of similar magnitude, 14 is in fair agreement 
with the measured H^ 2 (T). It is obvious that H p (T) docs 
not follow the same general T-dependence. Since thermal 
conductivity experiments probe the bulk of the sample, 
it is H C K {T) rather than H c p (T) that ought to be identi- 
fied as the upper critical field H^ 2 (T). Recent magneti- 
zation measurements 15 on single-crystals of MgB2 using 
a torque magnetometer result in values and a tempera- 
ture dependence of H^ 2 (T) that are consistent with our 
H%(T) and thus support our conclusion. Employing the 
equation given by anisotropic Ginzburg-Landau theory 
H c2 {9) = [(sm6/H^) 2 + (cos9/H^) 2 y 1/2 , where 9 is 
the angle between the magnetic field and the afr-plane, 16 
we estimate the errors in calculating H c2 caused by the 
above mentioned misalignment of 3.5 ± 0.5° to be about 
3 ± 1% for H^ and below 0.2% for H c c2 . 

As displayed in Fig. 4, at temperatures above about 



3 



H\\ab 



SAAaAa aaa Aa a 




31.3 K 



aaa a aaaa- 
■ n ou n oooo uuo 

28.2 K 



T = 28.2 K 



H increasing 
□ H decreasing 
O FC 



0.971— 
0.0 



0.90 



0.4 0.8 
H (kOe) 
I , L 



1.2 



10 



20 30 
H (kOe) 



40 




20 30 
T (K) 



40 



FIG. 4: Critical fields H p and H K = H C 2, as determined from 
electrical resistivity and thermal conductivity measurements, 
respectively. The solid line is compatible with calculations 
due to Helfand and Werthamer 14 . The dashed an dotted lines 
are to guide the eye. 



FIG. 3: The magnetic field dependence of the thermal con- 
ductivity n{H) for H _L c at T = 28.2 and 31.3 K. The solid 
vertical arrows mark H% b , set equal to the upper critical field 
Hc2- The inset demonstrates the small hysteresis of k(H) at 
low fields for T — 28.2 K. 



27 K, H^ 2 varies linearly with temperature with a 
slope dH^/dT = -1.17 kOe/K. This behavior leads to 
an extrapolated zero-field T' c = 36.6 K, 1.5 K below 
T c obtained from p(T, 0). Using the equations which 
are given by the Ginzburg-Landau theory considering 
anisotropies, 16 £ o6 (T) = ($ /2^ff c c 2 (T)) 1/2 and £ ab (T) = 
0.74(1 — T/T c )~ 1/,2 c; ah _ , where £ a (, is the coherence length 
in the basal a6-plane, we obtain the zero-temperature 
value of £ a fc.o = 11.8 nm. 

Turning to the temperature dependence of the critical 
field H^{T) for if _L c, we again note a sizable temper- 
ature interval where H^(T) varies linearly with T. This 
is emphasized by the dotted line in Fig. 4. The slope 
dH^/dT = -5.15 kOe/K gives V^otlo = 5.75 nm 
and therefore, the zero-temperature value of the c-axis 
correlation length £ Ci o = 2.8 nm. Another important pa- 
rameter which can be estimated from our k(H) data is 
the lower critical field if"*, as demonstrated in the in- 
set of Fig. 3. In the temperature region between 28 and 
35 K, where an evaluation of if?* with reasonable accu- 
racy of about ±10% was possible, H^/H^ « 130. From 
this ratio, using the equation if c2 /if c i — 2«q L / In kql, 16 
the parameter Kq L of the Ginzburg-Landau theory is es- 
timated to be about 13. 



As may be seen in Fig. 4, above approximately 33 K, 
if"! (f 1 ) deviates from the linear in T variation and, with 
increasing temperature, approaches zero also at T' c de- 
fined above. This is reflected in the temperature de- 
pendence of the anisotropy ratio 7, which seems to de- 
crease with T approaching T' c . The positive curvature 
of H C 2(T) is typical for strongly anisotropic, layered 
superconductors 17 and has often been explained in terms 
of the Lawrence-Doniach model, 18 which treats a lay- 
ered superconductor as a stacked array of weakly cou- 
pled two-dimensional superconducting sheets. Various 
other theoretical models have been proposed to explain 
this feature (for a critical review see, e.g., Ref. 19). At 
this point we cannot commit ourselves to any of these 
models. It is important, however, that the anomaly is 
absent for if || c and small and restricted to a rather 
narrow temperature region for if _L c. At lower tem- 
peratures, with decreasing temperature the anisotropy 
ratio 7(T) tends to a constant value and is approaching 
7o = H$(0)/H c c 2 (0) = &,6,o/£c,o = 4.2. 

The extrapolation to zero-temperature gives if^ 2 (0) w 
31 kOe, a considerably lower value than is typically 
claimed for MgB 2 . 5 Exceptions are the reports of 
Refs. 20,21. The rather low value of ff c c 2 (0) and the 
observation of a Helfand- Werthamer-type 14 temperature 
dependence of H^ 2 (T) have important consequences for 
possible models of the superconducting state of MgB 2 . 

Our result obviously questions the intrinsic nature of 
if C 2 (T) derived from measurements of p(T) . Since the re- 
sistive transition is not manifest in k(T), which may be 
considered as a bulk property, H p (T) must correspond 
to a minor fraction of an additional phase (or phases) 
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with enhanced H C 2 and T c . The spatial extension of 
this phase is, however, large enough to short-circuit the 
electrical current path and produce a narrow supercon- 
ducting transition at a temperature T c higher than the 
bulk transition temperature T' c . Based on an analysis of 
their magnetization and ac-susceptibility data on poly- 
crystalline samples, the authors of Ref. 22 came to a 
similar conclusion. The most likely origin of the second 
phase with enhanced superconducting parameters seems 
to be related to surface effects. A considerable enhance- 
ment of the electronic density of states near the Fermi 
level and, therefore, an enhanced trend to superconduc- 
tivity at the surface of MgB 2 have been predicted 23 ' 24,25 , 
in agreement with our observations. 

In conclusion, we observe a striking disagreement in 
the values and the temperature dependences of the up- 
per critical field H£ 2 of MgB 2 evaluated from results of 



electrical and thermal conductivity measurements on the 
same sample. The shape of H^ 2 (T) as established by 
k(H) with H || c does not reveal an anomalous posi- 
tive curvature near T c and therefore no exotic mechanism 
needs to be involved to explain the upper critical field, 
at least not for the bulk. Our data also indicate that the 
bulk transition temperature T' c is lower than T c obtained 
from results of p(T) measurements. 
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